Density Dependence of the Symmetry Energy and the Equation of State of Isospin 

Asymmetric Nuclear Matter 
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The density dependence of the symmetry energy in the equation of state of isospin asymmetric 
nuclear matter is studied using the isoscaling of the fragment yields and the antisymmetrized molec- 
ular dynamic calculation. It is observed that the experimental data at low densities are consistent 
with the form of symmetry energy, E syrn « 31.6 (p/po) ' 69 , in close agreement with those predicted 
by the results of variational many-body calculation. A comparison of the present result with those 
reported recently using the NSCL-MSU data suggests that the heavy ion studies favor a dependence 
of the form, E sym ~ 31.6 (p/p ) 7 , where 7 = 0.6 - 1.05. This constrains the form of the density 
dependence of the symmetry energy at higher densities, ruling out an extremely " stiff " and " soft 
" dependences. 

PACS numbers: 26.60.+C, 25.70.Pq, 25.70.Mn, 25.70.-z 



The Equation Of State (EOS) of isospin asymmetric 
(N Z) nuclear matter is a fundamental quantity that 
determines the properties of systems as small and light 
as an atomic nucleus, and as large and heavy as a neu- 
tron star P, H [1]. The key ingredient in the EOS of 
asymmetric nuclear matter is the density dependence of 
the symmetry energy. Theoretical studies [j, H H 0, [H 
based on microscopic many-body calculations and phe- 
nomenological approaches predict various different forms 
of the density dependence of the symmetry energy. In 
general, two different forms have been identified 0. One, 
where the symmetry energy increases monotonically with 
increasing density (" stiff " dependence) and the other, 
where the symmetry energy increases initially up to nor- 
mal nuclear density and then decreases at higher densities 
(" soft " dependence). 

Determining the exact form of the density dependence 
of the symmetry energy is important for studying the 
structure of neutron-rich nuclei [H, EH El, EH j and stud- 
ies relating to astrophysical origin, such as the structure 
of neutron stars and the dynamics of supernova collapse 
0, El El E3, El, El- For example, a " stiff " density 
dependence of the symmetry energy is predicted to lead 
to a large neutron skin thickness compared to a " soft " 
dependence p], El, HI El . Similarly, a " stiff" depen- 
dence of the symmetry energy can result in rapid cooling 
of a neutron star, and a larger neutron star radius, com- 
pared to a soft density dependence [H, [23| . 

In a heavy ion reaction, the dynamics of the collision 
between two heavy nuclei is also sensitive to the den- 
sity dependence of the symmetry energy [24], [25[. One 
can therefore carry out laboratory-based experiments to 
constrain this dependence. Recently 26], the fragment 
yields from heavy ion collisions simulated within the An- 
tisymmetrized Molecular Dynamics (AMD) calculation 
were reported to follow a scaling behavior of the type, 



Y 2 (N, Z)/Y 1 (N, Z) cx e aN+pz 



(1) 



are the yields from two different reactions. A linear rela- 
tion between the isoscaling parameter a, and the differ- 
ence in the isospin asymmetry (Z/A) 2 of the fragments, 
with appreciably different slopes, was predicted for two 
different forms of the density dependence of the symme- 
try energy ; a " stiff " dependence (obtained from Gogny- 
AS interaction) and a " soft " dependence (obtained from 
Gogny interaction). 

In this work, we show that the experimentally mea- 
sured scaling parameter a, favors a stiff density depen- 
dence of the symmetry energy, i.e. Gogny-AS interaction, 
and can be parametrized as E SJ(m w 31.6 (p/ Po) 1 , where 
7 = 0.69. The present observation is consistent with the 
EOS of Akmal and Pandharipande obtained from the 
many-body variational calculations [27], HI] . 

The measurements were carried out at the Cyclotron 
Institute, Texas A&M University using beams of 40 Ar, 
40 Ca, 58 Fe and 58 Ni from the K500 Superconducting Cy- 
clotron on 58 Fe and 58 Ni targets at 25, 30, 33, 40, 45, 47 
and 53 MeV/nucleon. Details of the experimental mea- 
surements and analysis can be found in Ref . [29J . 

Figure 1, shows the experimentally determined isoscal- 
ing parameter a, obtained from the fragment yields as a 
function of the beam energy. The different symbols cor- 
respond to various combinations of the reactions chosen 
for extracting the isoscaling parameters. The solid and 
the dotted lines are the exponential fits to the data. 

As mentioned earlier and shown in [26| , the parameter 
a, is related to the difference in the fragment isospin 
asymmetry (Z/A) 2 , through a linear relation of the form 



4C, 



T 



-[(Z/A) 2 -(Z/A) 2 



(2) 



where the parameters a and j3 are related to the neutron- 
proton content of the fragmenting source, and Y\ and Y<x 



where C sym is the symmetry energy and T is the tem- 
perature at which the fragments are formed. The quan- 
tity (Z/A) 2 — (Z/A) 2 ,, is the difference in the isospin 
asymmetry of the fragments in the two reaction systems. 
For the present systems, the isospin asymmetry of the 
fragments were evaluated at t = 300 fm/c of the dynam- 
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FIG. 1: Experimental isoscaling parameter a, as a function 
of the beam energy. The solid circles are from the Ar + Fc 
and Ca + Ni reactions. The open triangles are from Fe + Fc 
and Ni + Ni reactions. The solid stars are from Ar + Ni and 
Ca + Ni reactions. The open squares are from Fe + Ni and 
Ni + Ni reactions. 



ical evolution from the AMD calculations as discussed 
extensively in Rcf. [29| . 

Fig. 2 shows the a parameters plotted as a function of 
the difference in the fragment asymmetry for the beam 
energy of 35 MeV/nucleon. The solid and the dotted 
lines are the AMD predictions using the " soft " (Gogny) 
and the " stiff " (Gogny-AS) density dependence of the 
symmetry energy, respectively. The solid and the hollow 
symbols (squares, stars, triangles and circles) are the re- 
sults of the present measurements for the two different 
values of the fragment asymmetry, assuming Gogny and 
Gogny-AS interactions, respectively. Also shown in the 
figure are the scaling parameters (asterisks, crosses, di- 
amond and inverted triangle) taken from various other 
works [30l [3l| in the literature. One observes from the 
figure that the experimentally determined a parameter 
increases linearly with increasing difference in the asym- 
metry of the two systems as predicted by the AMD cal- 
culation. Also the data points are in closer agreement 
with those predicted by the Gogny-AS interaction (dot- 
ted line) than those from the usual Gogny force (solid 
line). 

It should be mentioned that in the above comparison 



FIG. 2: Scaling parameter a, as a function of the differ- 
ence in fragment asymmetry for 35 MeV/nucleon. The solid 
and the dotted lines are the AMD calculations for the Gogny 
and Gogny-AS interactions, respectively [2(|. The solid and 
the hollow, squares, stars, triangles and circles are from the 
present work as described in the text. The other symbols cor- 
responds to data taken from [3(j (asterisks) and [3lj ] (crosses, 
diamonds, inverted triangles). 



between the data and the calculation, the corrections 
for the isoscaling parameter a due to the secondary de- 
excitation of the fragments are not taken into account. 
The slightly lower values of the isoscaling parameters 
(symbols) from the present measurements with respect 
to the Gogny-AS values (dotted line) could be due to the 
small secondary de-excitation effect of the fragments not 
accounted for in this comparison. It has been reported 
by Ono et al. [32], that the sequential decay effect in 
the dynamical calculations can affect the a value by as 
much as 50 %. On the other hand, dynamical calculation 
carried out by Tian et al., [33| for the same systems and 
energy as studied by Ono et al., using Isospin Quantum 
Molecular Dynamic (IQMD) model, shown no significant 
difference between the primary and the secondary a. 

Due to the large discrepancy that exist in the deter- 
mination of the primary fragment excitation energy from 
the dynamical model calculations, it is difficult to esti- 
mate the effect of secondary de-excitation in dynamical 
models at this moment (34]. We have therefore assumed 
the effect of the sequential decay to be negligible, in the 
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above comparison. A small correction of about 10 - 15 
%, as determined from various statistical model studies 
[351 ] . results in a slight increase in the a values bringing 
them even closer to the dotted line. Note the asterisks 
symbols shown in the figure, and taken from the Ref. 
[301 ] . has already been corrected. The closer agreement 
of the experimental data with the Gogny-AS type of in- 
teraction, therefore, appears to suggest a stiffer density 
dependence of the symmetry energy rather than the soft 
Gogny interaction. 

Recently, Chen et al. [36] also showed, using 
the isospin dependent Boltzmann-Uehling-Uhlcnbcck 
(IBUU04) transport model calculation, that a stiff den- 
sity dependence of the symmetry energy parametrized 
as, E sym 31.6 (p/po) 1 ' 05 explains well the isospin dif- 
fusion data [S3| from NSCL-MSU (National Supercon- 
ducting Cyclotron Laboratory at Michigan State Univer- 
sity) . Their calculation was also based on a momentum- 
dependent Gogny effective interaction. However, the 
present measurements on isoscaling gives a slightly softer 
density dependence of the symmetry energy at higher 
densities than those obtained by Chen et al. 

This is clear from figure 3, which shows the parame- 
terization of various theoretical predictions of the density 
dependence of the nuclear symmetry energy in isospin 
asymmetric nuclear matter. The dot-dashed, dotted and 
the dashed curve corresponds to those from the momen- 
tum dependent Gogny interactions used by Chen et al. 
to explain the isospin diffusion data. These are given 
as, E S ym « 31.6 {p/poV, where, 7 = 1.6, 1.05 and 0.69, 
respectively. The solid curves and the solid points corre- 
sponds to those from the Gogny and Gogny-AS interac- 
tions used to compare with the present isoscaling data. 
As shown by Chen et al., the dependence parameterized 
by E sym w 31.6 [p/ p ) lm (dotted curve) explains the 
NSCL-MSU data on isospin diffusion quite well. On the 
other hand, the isoscaling data from the present work 
can be explained well by the Gogny-AS interaction (solid 
points). Both measurements yield similar results at low 
densities with significant difference at higher densities. It 
is interesting to note that by parameterizing the density 
dependence of the symmetry energy that explains the 
present isoscaling data, one gets, E SJ/m w 31.6 (p/po) J , 
where 7 = 0.69. This form of the density dependence 
of the symmetry energy is consistent with the parame- 
terization adopted by Heiselber g a nd Hjorth- Jensen in 
their studies on neutron stars [381 ] - By fitting earlier 
predictions of the variational calculations by Akmal et 
al. [27l . , where the many-body and special relativis- 
tic corrections are progressively incorporated, Heiselberg 
and Hjorth- Jensen obtained a value of E S j /m (p ) = 32 
MeV and 7 = 0.6, similar to those obtained from the 
present measurements. The present form of the den- 
sity dependence is also consistent with the findings of 
Khoa et al. [39j . where a comparison of the experimen- 
tal cross-sections in a charge-exchange reaction with the 
Hartree-Fock calculation using the CDM3Y6 interaction 
[40| . reproduces well the empirical half-density point of 
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FIG. 3: Parameterization of the various forms of the den- 
sity dependence of the nuclear symmetry energy used in the 
analysis of the present measurements on isoscaling data and 
the isospin diffusion measurements of NSCL-MSU ;36j. The 
various curves are as described in the text. 



the symmetry energy obtained from the present work (see 
fig. 2 of Ref. [l^). 

The observed difference in the form of the density de- 
pendence of the symmetry energy between the present 
measurement and those obtained by Chen et al. is not 
surprising. Both measurements probe the low density 
part of the symmetry energy and are thus less sensitive 
to the high density region. But the important point to be 
noted is that both measurements clearly favor a stiff den- 
sity dependence of the symmetry energy at higher densi- 
ties, ruling out the very " stiff " (dot-dashed curve) and 
very " soft " (solid curve) predictions. These results can 
thus be used to constrain the form of the density depen- 
dence of the symmetry energy at supranormal densities 
relevant for the neutron star studies. 

It should be mentioned that the calculations in both 
the above described works assume a similar value for 
the symmetry energy at normal nuclear density (about 
31 MeV) . Although, numerous many-body calculations 
d, [U [42|, 51] and those from the empirical liquid drop 
mass formula pit |45| predict symmetry energy near nor- 
mal nuclear density to be around 30 MeV, a direct ex- 
perimental determination of the symmetry energy does 
not exist. 
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Recently, Khoa et al. 13911 . analyzed the experimental 
cross-section data 0, |47| using the isospin dependent 
CDM3Y6 interaction of the optical potential in a charge 
exchange p( 6 i?e, 6 Li*)n reaction. Their analysis probed 
mainly the surface part of the form factor and hence ap- 
propriate for densities close to the normal nuclear density. 
Based on their results and the Hartree-Fock calculation 
of asymmetric nuclear matter using the same effective 
nucleon-nucleon interaction, they estimate the most re- 
alistic value of the symmetry energy to be about 31 MeV. 
An accurate determination of the neutron skin thickness 
AR, from the parity-violating electron scattering mea- 
surement 20] is however, likely to provide a more precise 
determination of the symmetry energy near normal nu- 
clear density. 

In view of the findings from the present measurements 
and those of Chen et al., we believe that the best esti- 
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